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We study level correlations in a two-dimensional system with a long-range random potential and 
strong spin-orbit (SO) splitting Aso of the spectrum. The level correlations for sufficiently large Aso 
are shown to be described by orthogonal Wigner-Dyson statistics, in contrast to the common point 
of view that a system with the strong SO coupling belongs to the symplectic statistical ensemble. 
We demonstrate also that in a wide energy interval the level statistics is completely determined 
by transitions between two branches of the split spectrum. A sharp resonance is obtained in the 
two-level correlation function at energy equal to Aso. 

PACS numbers: 73.20 Fz, 71.25 Mg, 72.15 Rn 
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The advancement of fabrication of nanostructure de- 
vices in the last two decades has initiated great interest 
to properties of small mcsoscopic systems (see Ref. ||] 
for review) . One of the most frequently studied charac- 
teristics of the mesoscopic samples is the two-level cor- 
relation function i?(e) = {p{E -\- e)p{E))/p'^ — 1 which 
describes fluctuations of the density of states p{E). Here 
p = {p{E)) and (...) denotes averaging over realizations of 
the random potential. The problem of level correlations 
in small disordered conductors was first discussed in pj, 
where it was suggested that the distribution of energy 
levels is determined by one of the three types of the uni- 
versal Wigner-Dyson statistics O] . Later on this idea was 
confirmed by rigorous calculations |J] of i?(e). It has been 
proved that the universal Wigner-Dyson distribution is 
realized at small energies e <^ E,., where E^ ^ 'hJT^ is 
the Thouless energy (here t^ ~ 1? jD, L is the size of 
the system and D is the diffusion coefficient). At larger 
energies (e ^ E^) the level statistics is nonuniversal and 
depends on the dimensionality of the system g. The 
most interesting is the 2d case in which the calculation 
of -R(e) in the standard one-loop approximation gives a 
zero result [|| . In this situation the level statistics is gov- 
erned by weak localization effects [Q and is sensitive to 
the boundary conditions and topology of the surface . 

In this paper we study level correlations in a 2d sys- 
tem with strong SO spectrum splitting A^-o, which orig- 
inates from an asymmetry of the confining potential of 
the quantum well [g[. The enhanced interest to the sys- 
tems with strong SO coupling was initiated by recent 
remarkable experiments of Kravchenko et al. [p| which 
clearly demonstrated the existence of a metal-insulator 
transition in 2d systems. The importance of the SO spec- 
trum splitting for the transition was first conjectured in 
Ref. pO| . This idea was strongly supported by very re- 
cent experimental data |0,fi3]. 

Here, we show that strong SO spectrum splitting leads 
to unexpected effects in the statistics of levels in a 2d dis- 
ordered system. A common belief is that the type of the 
level statistics in a system with the SO coupling depends 



only on the relation between the spin relaxation time Tso 
and Tc. For the case of very small splitting, when t^q > r^, 
spin degrees of freedom do not play any role and the sys- 
tem belongs to the orthogonal ensemble - one of the three 
universal Wigner-Dyson ensembles. With increasing SO 
coupling, Tso becomes smaller than Tc^ the spin relax- 
ation comes into play and the system exhibits crossover 
to symplectic ensemble. This scenario takes place for the 
case of a short-range potential. As we show below, the 
system with long range disorder demonstrates reentrance 
to the orthogonal ensemble at sufficiently large values of 
Aso- This surprising result is related to strong suppres- 
sion of the transitions between two branches of the SO 
split spectrum in the case Aso ^ fivp/d [^ (here d is 
the potential correlation length, vp is the Fermi velocity). 
As a consequence, the system is divided into two inde- 
pendent subsystems, the spin degrees of freedom are ef- 
fectively frozen out, so that in the limit of very large Ago 
we have two independent orthogonal ensembles. We will 
demonstrate also that long-range nature of the impurity 
potential is of crucial importance for the level statistics 
in the nonuniversal region of energies e ^ Ec- We argue 
that in a wide energy interval in the nonuniversal region 
the level correlations are completely governed by the rate 
of transition between two branches of the spectrum and 
that a narrow resonant peak arises in R{e) for e ~ Ago- 

The term in the Hamiltonian which is responsible for 
the SO splitting of the energy spectrum B is given by 
a{d- X k) • n, where a are the Pauli matrices, Jik is the 
electron momentum and n is the unit vector normal to 
the 2d layer. A constant a characterizing the strength 
of the SO coupling is proportional to the electric field 
arising due to asymmetry of the confining potential of 
the quantum well. The value of a can be tuned by 
applying an external electric field perpendicular to the 
2d plane uM- The energy spectrum and spin wave- 
functions of an electron with the effective mass m read: 
E^{k) = h^P/2m±ak, Xk = (l,±ie*'^'') /V2. The 
spinors Xu. depend on the polar angle (pk of the wave vec- 
tor k and describe two states with spins polarized along 



the vectors ±(k x n), respectively. Thus the system is 
divided into two branches, (+) and (— ) (see Fig. la) and 
the electron spin in each of the branches is rigidly con- 
nected to the momentum. 

The long-range random impurity potential U{r) is as- 
sumed to be weak enough, so that the inequality r » 
d/vp holds, where t ~ Ttr/ikpdY is the elastic scat- 
tering time, kp is the Fermi wave- vector. The presence 
of this potential leads to transitions both within each 
branch and between different branches. The respective 
(intrabranch and interbranch) times of these transitions 
are given by the following expression 
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J^,,.(k-k')5[i?''(k)-i?''(k')], (1) 



where X^^(k-k') = |(x(;|xS:,)P^(|k-k'|), indices ^i and 
V denote the type of the branch (+ or — ), and the func- 
tion K{k) is the Fourier-transform of the potential corre- 
lation function K{r) = (C/(r)[/(0)). The function K{r) 
falls off on the scale of d. For the case when the spin-orbit 
splitting A^o = 2akp is not too large (Ago <C Ep, where 
Ep is the Fermi energy), we can set t+_|_ = r__ — t, 

T^ = T y- = T*. The minimal transferred momentum 

needed for the transition between the branches is equal 
to ImajTi. This momentum should be compared with 
?i/d. If the inequality A^o ^ Tivp /d is fulfilled, then 
the interbranch transitions are suppressed compared to 
transitions within one branch (note that this inequality 
implies also that A^oT ^ K). In particular, for the case 
of the potential created by ionized impurities located at 
distance d from the 2d layer, the correlation function 
K{q) ~ exp(— 2g(i) and the interbranch transition time 
is given by 



r, = 4:T{kpd)'^ exp 
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The factor (kpd)'^ in this expression is due to the orthog- 
onality of the spinors corresponding to different branches 
with the same direction of momenta. Note that t^ ~ Ttr 
for the case of a short-range potential {K{q) — const). 

Expressing the density of states in the usual 
way in terms of exact Green functions p{E) = 
{i/2TrL'^)TT{G"{E) - G^{E)), we can write the two-level 
correlation function as i?(e) = (AV87r2)Re(TrG^(^ -|- 
e)TrG'^ (£■)), where A = 27rh'^/mL^ is the mean level 
spacing and the trace is calculated over spatial and spin 
variables. The main contribution to the function R{f.) 
comes from one-loop diagrams containing two diffuson 
or cooperon propagators, shown in Fig. 2 in the so-called 
Hikami representation |lj|. In what follows we consider 
the case of zero magnetic field, so that the cooperon con- 
tribution is equal to the diffuson one. This allows us to 
focus on the calculation of the diffuson contribution only. 
To separate the two branches of the spectrum, we use the 
following representation for the averaged Green functions 



(G«.^(k, E)) = E^^± G«'^(k, E)P^{-k), where P^(k) = 
|x(^)(x(^| is the projection operator to the branch ^. The 
functions G^'"^ are represented by edges of the "Hikami 
boxes" in Fig. 2. It is convenient to distinguish two 
channels: the symmetric channel with fj, = v, ^' — v' 
and antisymmetric channel with /i = — i^, p! = —v' (see 
Fig. 2). We will first concentrate on the region of small 
energies {tTtr *C ?i), where the symmetric contribution 
dominates. The diffuson propagator in symmetric chan- 
nel depends on two indices ^ and ^' only: D^^, = D^^ . 
This propagator corresponds to the contribution of two 
coherent wave ("particle" and "hole" waves) which be- 
long to the same branch (see Fig. lb) in any part of 
the diffusive trajectory. The value Z?^'^ is determined 
from the system of four equations (/i — ±,/j,' = ±) 
depicted graphically in Fig. 3, where X^^ (k, k',q) = 

if (|k-k'|)(xi:ixi:;) (xC;+qlxC+q)- These equations should 
be solved together with the system of equations for 
G^'^(k, i?) written in the framework of Self-Consistent 
Born Approximation. In the diffusion approximation, the 
solution is given by D++(q, e) — D (q, e) = S'+Ti/Trpr^ 
and L>+-(q,e) ^ D-+{ci,e) ^ S-h/npT^ ||, where 
functions 
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obey the set of equations fl^ 
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dt 



DAS+ ^ S{r)S{t) + 
DAS_ = ^+-^- . 
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Here D = vpTtr/2 is the diffusion coefficient, F* = 2h/T*. 
Note that in the case of a short-range potential r* ^^ Ttr 
and the term {HDq'^ — ie + T^,)^^ can be neglected in both 
5"+ and S- . The contribution of the symmetric channel 
to the two-level correlation function is given by 
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The first term in the r.h.s of this expression is the well- 
known Altshuler-Shklovskii result H. The second term 
is related to the interbranch transitions. 

First we discuss the universal region of energies, e <C 
Ec- For r* ^ Ec Pql one should keep only the term 
q = in Eq. (pi) , which yields 
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A crossover between two ensembles is clearly seen. For 
e <C r* the first term in Eq. m) dominates. In this 



case during the time U/e many interbranch transitions 
take place, the system is ergodic and we obtain the stan- 
dard result [H for the envelope of the symplectic Wigner- 
Dyson distribution. For F* <C e <C i^^c the time ~ h/t is 
not enough for interbranch transition and we have two 
independent orthogonal ensembles. Thus, we see that the 
system demonstrates reentrance to the orthogonal ensem- 
ble at large A^o, when F* becomes sufficiently small. 

In the nonuniversal region e '^ Ec, where the summa- 
tion in Eq. (0) can be replaced by the integration, the 
contribution of the first term in Eq. (H) is equal to zero 
(for periodic boundary conditions) and the level correla- 
tions are governed by the interbranch transitions only: 






(7) 



Here g = irkpltr is the conductance in units of c^/tt/i. 

Compare now Eq. (R) with the weak localization cor- 
rection to i?(e) for the symplectic ensemble [^ Rwi{^) = 
A/Sg'^e. The inequality Rs{e) > Rn,i{e) is fulfilled in the 
energy interval max(_Ec,F*/(7) <^ e <C g^*- It can be 
shown |l7| that the boundary-induced and topological 
contributions H are also small compared with Eq. ffl). 
Thus, we see that in the nonuniversal region the inter- 
branch transitions completely govern level correlation in 
a wide energy interval. 

Now we will show that in the region of large energies 
\e\ ~ Aso a narrow resonant peak appears in the two- level 
correlation function. The physical origin of this peak is 
this: the exact energy levels lying near the energy E in 
the lower branch (— ) are correlated with the exact en- 
ergy levels lying near the energy E + Ago in the upper 
branch (-I-) (see Fig. Ic). The resonance contribution to 
R{e) comes from the diagram in Fig. 2 with /i = —v 
and /i' = —v' (antisymmetric channel). The diffusons 
in antisymmetric channel depends on indices fi and n' : 
D!!/^L^, = D^^''. The equations for D^^ are obtained 
from equations for D^^ by changing the sign of all the 
indices related to the line corresponding G (lower line 
in Fig. 3). The dashed line represents now the func- 

tion Kr'(k,k',q) = K{\k~k'\){x^\xt,){Xu''UxZ.,)- 
In the case e « A^o the main contribution is represented 
byl?++(q,e) = (h^ /npT^)/ {hDq^ + T-i[e- Aso]), where 
F — Ti/t and 
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Here Ktrik - k') = /Xdk - k'|)(l - cos{ip))dLp and ^ 
is the angle between k and k'. At large spectrum split- 
ting the inequalities r^^ <C t ^ r* hold. Now we present 
qualitative arguments clarifying the physical origin of r 
which plays the role of the phase breaking time in the 
antisymmetric channel. The appearance of a resonance 
peak in i?(e) is a consequence of coherence of the waves 



(-I-) and (— ), which for e = A^o have the same momenta 
(see Fig. Ic). One might think that the coherence is de- 
stroyed by interbranch transitions only. However, in the 
second order of perturbation theory there exists a vir- 
tual process which breaks the coherence. In this process, 
the wave which in the initial and final state belongs to 
the branch (+), in the intermediate virtual state belongs 
to the branch (— ). This process leads to a small cor- 
rection to the effective random potential which can be 
estimated as VF ~ U"^ / AsoikpdY (here U is the ampli- 
tude of random potential and the factor {kpd)~'^ is of 
the same origin as in Eq. (ph). The other wave, in the 
initial and final state belonging to the branch (— ), in the 
intermediate virtual state belongs to the branch (-I-). It 
is evident that the correction to the potential W has the 
same absolute value but the different sign W' = —W. 
The sign difference is related to the fact that the energy 
of the virtual state lies in first process above, and in the 
second process below energy of the initial state. As a 
result, two waves propagate in slightly different effective 
potentials U + W and U ~ W. On passing the poten- 
tial correlation length d they acquire the phase difference 
S(j) r^ Wd/vpTi- The phase breaking time is then easily 
estimated as f~^ ~ {Scfy^vp/d. Combining these formu- 
las one can obtain the estimate for f given in the r.h.s. 
of Eq. (0). The resonance contribution to R{e) reads 
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(9) 



Sir' [(e - A,o)2 + F2 



A < F < ^e 



We conclude that a 2d system with a long-range ran- 
dom potential and strong SO splitting of the spectrum 
exhibits a reentrance to the orthogonal ensemble for suf- 
ficiently large values of A^^. A weak interaction between 
two sets of random levels corresponding to two branches 
of the split spectrum leads to the appearance of two ef- 
fects: (i) in a wide energy interval the level statistics 
is completely governed by interbranch transitions; (ii) a 
sharp resonance arises in the two-level correlation func- 
tion for energies corresponding to the spectrum splitting. 
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FIGURE CAPTIONS 

Fig. 1 (a) Spin-orbit split energy spectrum; (b) sym- 
metric contribution to the difFuson (two coherent waves 
are depicted by black dots); (c) antisymmetric contribu- 
tions to the difFuson. 

Fig. 2 One- loop diagram contributing to -R(e). 

Fig. 3 DifFuson equation in the symmetric channel. 
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